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Abstract 

We analyze the production of neutral Higgs particles in association with the su- 
persymmetric scalar partners of the third generation quarks at future high-energy 
hadron colliders [upgraded Tevatron, LHC] and e + e~ linear machines [including 
the 77 option]. In the Minimal Supersymmetric extension of the Standard Model, 
the cross section for the associated production of the lightest neutral h boson with 
the lightest top squark pairs can be rather substantial at high energies. This pro- 
cess would open a window for the measurement of the htt coupling, the potentially 
largest coupling in the supersymmetric theory. 
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1. Introduction 



One of the firm predictions of the Minimal Supersymmetric Standard Model (MSSM) Q 
is that, among the quintet of scalar [two CP-even h and H, a pseudoscalar A and two 
charged if ± ] particles contained in the extended Higgs sector |2j , the lightest Higgs boson 
h should be rather light, with a mass Mh < 130 GeV [||, |J. This Higgs particle should 
therefore be produced at the next generation of future high-energy colliders ||, Q , or even 
before [Tj, |8[], if supersymmetry (at least its minimal version) is indeed realized in Nature. 
At the Large Hadron Collider (LHC), the most promising channel [[| for detecting the h 
boson is the rare decay into two photons, with the Higgs particle dominantly produced via 
the gluon fusion mechanism ||; the two LHC collaborations expect to detect the narrow 
77 peak for M h < 130 GeV, with an integrated luminosity / Cdt ~ 300 fb _1 0. At future 
high-energy and high-luminosity e + e~ colliders, with center of mass energies in the range 
350-500 GeV, the properties of this particle should be studied in great detail ||, to shed 
some light on the electroweak symmetry breaking mechanism. 

In the MSSM, the genuine SUSY particles, neutralino/charginos and sfermions could 
have masses not too far from the electroweak symmetry breaking scale. In particular 
the lightest neutralino, which is expected to be the lightest SUSY particle (LSP), could 
have a mass in the range of ~ 100 GeV. Another particle which could also be light is 
one of the spin-zero partners of the top quark, the lightest stop t±. Indeed, because of 
the large m t value, the two current stop eigenstates could strongly mix [[l(J, leading to 



a mass eigenstate t\ much lighter than the other squarks [which are constrained to be 
rather heavy [jTTJ by the negative searches at the Tevatron] and even lighter than the top 
quark itself. Similar features can also occur in the sbottom sector. These particles could 
therefore be also easily accessible at the next generation of hadron and e + e~ colliders. 

If the mixing between third generation squarks is large, stops/sbottoms can not only 
be rather light but at the same time their couplings to Higgs bosons can be strongly 
enhanced. In particular, the (normalized) htit\ coupling can be the largest electroweak 
coupling in the MSSM. This might have a rather large impact on the phenomenology of 
the MSSM Higgs bosons as was stressed in Ref. [fT5| . The measurement of this important 



coupling would open a window to probe directly some of the soft-SUSY breaking terms 
of the potential. To measure Higgs-squarks couplings directly, one needs to consider the 
three-body associated production of Higgs bosons with scalar quark pairs which has been 
studied recently in [TJ|, [L5L This is similar to the processes of Higgs boson radiation 



from the top quark lines at hadron [[17]] and e + e [18[ colliders, which allow to probe the 
tt-Higgs Yukawa coupling directly. 



In this paper we extend our previous study |14j and discuss the production of neutral 



Higgs particles in association with third generation scalar quarks at future high-energy 
hadron and e + e~ linear machines [including the 77 option] in the unconstrained MSSM 
and minimal SUGRA. In the next section, we will summarize the properties of third 
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generation squarks and their couplings [including the constraints on the latter]. In section 
3 and 4, we discuss the associated production of squarks and Higgs bosons at hadron and 
e + e~ colliders, respectively. Section 5 is devoted to the conclusions. 

2. The Physical Set-Up 

There are two main reasons which make the rate for the associated production of Higgs 
bosons with third generation scalar quarks potentially substantial: due to the large 
Yukawa couplings, the lightest top and bottom squarks can have relatively small masses, 
and their couplings to Higgs bosons can possibly be large. In this section, we will first 
summarize the properties of top and bottom squarks and discuss their masses, mixing and 
their couplings to the MSSM neutral Higgs bosons [and to the massive gauge bosons], and 
then list the various experimental and theoretical constraints which will be taken into ac- 
count in our analysis. 

2.1 Masses and couplings 

In the case of the third generation, the left- and right-handed sfermions q~L and q R [the 
current eigenstates] can strongly mix |T(J ; for a given squark q = t,b, the mass matrices 



which determine the mixing are given by 

<m ^ nm A 

(1) 



rn 2 T r, + m 2 m g A 



Mh m 2 RR + m 2 

with, in terms of the soft SUSY-breaking scalar masses m,q L and m^, the trilinear squark 
coupling A q , the higgsino mass parameter /i and tan/? = vu/vd, the ratio of the vacuum 
expectation values of the two-Higgs doublet fields: 

m LL = m \ L + ( J 3 - e q s w) COS W M Z , m RR = m \ R + e q s W COS 2/5 M Z 

A q = Aq- f i(t ( mf3)- 2I $ (2) 

If and e q are the weak isospin and electric charge (in units of the electron charge) of the 
squark q, and s^y — 1 — = sin 2 #iy- The mass matrices are diagonalized by 2 x 2 
rotation matrices of angle 9 q 

B? = f Cdq Sdq | , Co = cos a and sg„ = sin a (3) 
V s dq c 0q J 

The mixing angles 9 q and the squark eigenstate masses are then given by 

2m A a „ m 2 T r — m% R , . 

sin 20, = 9 q S > cos2d g = — F 4 

m\ x - m| 2 ml - m\ 2 



i i 
m- = m„ H — 
91,2 q 2 



m LL + m RR =F V ( m LL ~ ™RR) 2 + ^\A\ 



(5) 
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Due to the large value of m t , the mixing is particularly strong in the stop sector, unless 
A t = 0. This generates a large splitting between the masses of the two stop eigenstates, 
possibly leading to a lightest top squark much lighter than the other squarks and even 
lighter than the top quark. For large values of tan/3, /x and/or A\, the mixing in the 
sbottom sector can also be rather large, leading in this case to a possibly light b\. 



In the constrained MSSM or minimal SUGRA [19], the soft SUSY breaking scalar 
masses, gaugino masses and trilinear couplings are universal at the GUT scale; the left— 
and right-handed sfermion masses are then given in terms of the gaugino mass parameter 
mi/2, the universal scalar mass m and tan/3. For instance in the small tan/3 regime, 
the two soft top scalar masses [in which we include also the D-terms] at the low energy 
scale obtained from the one-loop renormalization group evolution ||20|| , are approximately 
given by 



m\ L = 0.51m^ + 5.29m2 /2 + 0.35M|cos2/3 

m\ R = 0.01m2 + 3.68m2 /2 + 0.15M|cos2/3 (6) 

This shows that, in contrast with the first two generations, one has generically a sizeable 
splitting between m~ L and at the electroweak scale, due to the running of the (large) 
top Yukawa coupling [expect possibly in the marginal region where m and both 
become very small]. Thus degenerate left and right soft stop masses [at the electroweak 
scale] will be considered only in the unconstrained MSSM. Note that in this case, in most 
of the MSSM parameter space, the mixing angle eq. @ is either close to zero [no mixing] 
or to — 7r/4 [maximal mixing] for respectively, small and large values of the off-diagonal 
entry m t A t of the squark mass matrix eq. (1). 

The couplings of a neutral MSSM Higgs boson $ = h, H, A to a pair of squarks 
is given by 

= £ ca, (*% (7) 

k,l=l J 

with the matrices Cq>qq summarizing the couplings of the Higgs bosons to the squark 
current eigenstates; normalized to 2(v2G f) 1 ^ 2 Mz they are given by 

Ch __ , " (I! - e q * 2 w) sin(/3 + a) + %s\ ^{A q s\ + ^) ^ 
7^r(A q s\ + ns q 2 ) -e q s 2 w sin(/3 + a) + 



C Hm = \ m„ , , a . a, , , n . s . ml „ \ ( 9 ) 



q_ r 1 

1 1 
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CAqq 



1 

2 



Ml 



Mi 



H + A, (tan [3)- 2L * 



(i + A g (tan/3)- 2J 3 




(10) 



with a the mixing angle in the CP-even Higgs sector, and the coefficients r\ 2 and s\ 2 for 
top and bottom squarks given by 



' 2 



COS OL 



4 



sin a 



sin a 



4 



cos a 



11^ 



sin/5 ' " 1 sin/3 ' 1 2 cos/5 ' 2 1 cos/? 

In the decoupling limit, Ma ~ M# 3> M z , the mixing angle a reaches the limit (3 — n/2, 
and the expressions of the couplings become rather simple. One would then have for the 
hit couplings for instance 



9 huh 

9ht 2 t 2 



cos 2(3 
cos 2(3 
cos 2(3 sin 29 t 



cos 2 6 t — 



sin 2 6 t 



s w cos 26> £ 



+ 



m 



M 2 



+ - sin 26 t 

z 2 



s w cos 26 1 



+ 



m 



M 2 



m t A t 
M 2 
m t A t 



Sin2 ^ Ml 



1 m t i4 ( 
+ 2 COS2 ^MF 



(12) 



and involve components which are proportional to A t — A t — /i/ tan/3. For large values of 
the parameter A t which incidentally make the t mixing angle almost maximal, | sin2^| ~ 
1, the last components can strongly enhance the coupling and make it larger than 

the top quark coupling of the h boson, ghtt oc m t /Mz- The couplings of the heavy H 
boson to stops involve also components which can be large; in the case of the lightest 
stop, the coupling reads in the decoupling limit: 

,2 i 



cos &i — -s w cos 26 1 



m; 



1 „ TUf 

- sin 26 
2 



At 



Mrtan/3 



+ /i) (13) 



Ml tun (3 & iv± z 

For large tan (3 values, the m 2 and the A t components of the Qhi^ coupling are suppressed; 
only the component proportional to \i is untouched. The pseudoscalar A boson couples 
only to ti?2 pairs because of CP-invariance, the coupling is given by: 

1 m t 



5 At lt 2 



2 Ml 



(A t / tan P + fi) 



(14) 



In the maximal mixing case, | sin2# t ~| ~ 1, this is also the main component of the H boson 
coupling to tit 2 pairs, except that the sign of \i is reversed. 

Finally, we will need in our analysis the squark couplings to the Z boson. In the 
presence of mixing and normalized to the electric charge, the coupling V z i i i j = eaij, are 
given by 



2-^3C# 2Su/Ca 2I%S 2 2Su/6 n 



a u = 



2c\yS\y 



a 2 2 = 



2C\\rS\\r 
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a 12 = a 2 i = 



2c\yS\y 



(15) 



The couplings to the photons, with the same normalization, are simply given by eqSij. 
With this notation the vector and axial-vector couplings of the electron to the Z boson, 
that we will also need later, read [with the vertex convention e^^{y e — a e 7 5 )] 

-_-!+<+ i (16) 



2.2 Constraints used in the numerical analysis 

We will illustrate our numerical results in the case of the unconstrained MSSM as well as 
in the minimal SUGRA case. We will concentrate on the case of the lightest h boson of the 
MSSM in the decoupling regime [^| , and discuss only the production in association with 
light top squarks [in the case of no mixing, this corresponds to both squarks ti,^ since 
their masses are almost equal, while in the case of large mixing this would correspond to 
the lightest top squark ti\. That is, we consider the range of the MSSM input parameters 
in such a way that Mh <C Mh ~ Ma ~ M#±, which also implies for the mixing angle a 
of the neutral CP-even Higgs sector a ~ (3 — tt/2. In this approximation, the behavior 
of the tth coupling with respect to the SUSY-breaking parameter A t is very simple, as 
discussed in the previous subsection. 

The approximation of being close to the decoupling limit implies that we do not 
consider other Higgs production processes [although our analytical expressions will hold 
in these cases] with e.g. the heavy neutral H or CP-odd A bosons produced in association 
with top squarks. In any case, these processes would be suppressed by phase-space well 
before the decoupling regime is reached [this is even more legitimate in the minimal 
SUGRA case, where for most of the parameter space allowed by present experimental 
limits, the MSSM Higgs sector turns out to be in this decoupling regime; see Ref. []23| 
for instance.] The production of the h boson in association with light sbottoms [that we 
will also not consider in the numerical analysis] should follow the same pattern as for the 
stops, except that one has to replace the t\t\h coupling by the bibih coupling; the latter 
being strongly enhanced for large tan j3 and \x values. 

Note that, in general, we will use as input the parameters rrig,A t , tan/3 and \x. The 
h boson mass is then calculated as a function of these parameters [Mh is only marginally 
affected by the variation of the parameter u] with the pseudoscalar Higgs boson mass 
fixed to Ma = 1 TeV, and with the full radiative corrections in the improved effective 
potential approach || included. 

Let us now discuss the constraints which will be used in our numerical analysis. First 
there are of course the experimental bounds on the top squark and Higgs boson masses. 
From negative searches at LEP2, a bound > 80 GeV is set on the lightest top squark 
11]. Larger masses, up to m~ tl ~ 120 GeV, are excluded by the CDF collaboration 



mass 



Eg] in the case where the lightest neutralino is not too heavy, m x o < 50 GeV. 
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The general experimental lower bound on the lightest h boson in the MSSM from 
LEP2 negative searches is > 80 GeV jO]]. However, in the decoupling limit the h 



boson will behave as the SM Higgs particle for which the lower bound is larger, > 90 
GeV; this will be the lower bound which will be used practically in our analysis since we 
will work mostly in the decoupling regime. This calls for relatively large values of the 
parameter tan/3 [since the maximal value of Mh increases with increasing tan/5] and/or 
large mixing in the top squark sector [the "maximal mixing" scenario, in which the h 
boson mass is maximized for a given tan (3 value, is for A t ~ y/Erriq]. 

In this paper, we will often focus on the large At/ 'small scenario, which gives 
the largest t±tih production cross sections. However there are theoretical as well as ex- 
perimental constraints on this scenario which should also be taken into account in the 
analysis: 

i) The absence of charge and color breaking minima (CCB) p4j can put rather strin- 
gent bounds on the parameter A t : for the unconstrained MSSM stringent CCB 



constraint for a large A t scenario reads |24 



A\ < 3{ml + m\ R + /i 2 + m% v ) (17) 
to be valid at the electroweak scale. [Here we only aim at a qualitative treatment and 



do not enter into possible improvements of condition (|T7j), see the last reference in [24 



Note, however, that for the large A t values we consider, the CCB condition from the 



precise comparison of the actual minima is practically equivalent to ( |TTD .] This constraint 
involves the stop sector parameters as defined in the previous subsection, plus m% u , the 
mass squared term of the scalar doublet $fy. Since, as indicated above, we shall fix in 
our analysis /i and tan/3 (and Ma — 1 TeV for the decoupling regime), m| is not an 
independent parameter and will be determined consistently from the electroweak breaking 
constraints. It can be easily seen that the constraint (pj]) is in fact weakly dependent on 
\i [and we also note as a general tendency that ( |TTD is slightly less restrictive for relatively 
small tan/5 values, tan/5 ~2-5, being however almost independent of tan/5 for tan/3 > 5]. 
Clearly thus, large values of A t call for larger values of m^, m i R and therefore generally 
larger values of mi x . For instance, in the mSUGRA case with the universality condition 
mi = mi = mo at the GUT scale, A t should be restricted to values smaller than 
~ 1.5(1.2) TeV for tan/3 ~ 3(30) and m h ~ 200 GeV. In the unconstrained MSSM case, 
one may take relatively large A t values by pushing e.g. m~ tR to larger values, while still 
allowing a relatively light tx, depending on the values of /i and tan/3. The price to pay, 
however, is that the mi -mi splitting increases, therefore disfavoring too large a value 



for the t]t-Ji coupling, as is clear from eqs. (£|) and ([12]) . Eq. (|17|) thus serves as a relevant 
"garde-fou" to anticipate in which region of the parameter space we can expect the qq§ 
cross sections to be substantial without conflicting with the CCB constraints. In our 
numerical illustrations in the next section, we shall indicate more precisely whenever the 
CCB constraint (|T7|) becomes restrictive. We note however that, after having delineated 
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the most dangerous CCB regions, quite often the CCB constraints are superseded either 
by the lower bound on the h boson mass and/or by the present indirect experimental 
constraints coming from the virtual contributions of squarks to the electroweak observables 
at the Z peak. Indeed, in general the latter also forbid very large values of the parameter 
At, as we shall discuss in more detail now. 

ii) Large values of A t lead to a large splitting of the top squark masses; this breaks 
the custodial SU(2) symmetry, generating potentially large contributions to electroweak 
high-precision observables, and in particular to the p parameter [^5[ , whose leading con- 
tributions are proportional to different combinations of the squared squark masses, and 



which is severely constrained by LEP1 data | 2"6|| . Since this is potentially the most restric- 



tive constraint on a large A t scenario, we have thus evaluated precisely the contributions 
of the (t, b) isodoublet to Ap = p — 1 for typical choices of the parameters corresponding 
to our illustrative examples [e.g. to Figs. 2 and 3 of next section]. Explicitly p5f : 



Ap(t,b) 



Gi 



•57T 



V2 



[cl (eg f(ml,ml) + si /(m|, m|)) + s 2 t (c 2 b /(m|, m£) + s 2 b /(m|, mj,)) 



-c 



\s\ f(m{,ml 2 ) - cist (18) 

where f(x,y) = x + y — 2xy/(x — y) ln(x/y) and s i: Ci designate the sine and cosine of 
the mixing angles in the stop and sbottom sectors. As expected, these contributions can 
become unacceptably large for very large A t and small m^, but stay however below the 
acceptable level, Ap\^- b < 3 ■ 1CT 3 [which approximately corresponds to a 2a deviation 
from the SM expectation ] for values of A t as large as 1.4-1.5 TeV and stop masses 
as small as > 130-140 GeV [those precise limit values depend also on the values 
of tan (3 and p] . We refrained however from making a systematic analysis of these low- 
energy constraints, since these effects are also quite dependent on the other parameters 
and in our numerical analysis we shall simply indicate whenever the contributions to the 
p parameter are becoming restrictive. We note, however, that the magnitude of Ap(t,b) 
in eq. ([18]) is essentially driven by the relative values of m t At and mq L : Ap(t,b) can 



be large when m t A t ~ m Q L while there is a decoupling, Ap(t,b) ~ l/m^, whenever 
m QL m tAt- (This decoupling is expected, since for very large vtlq l there is no more 
difference between the W± and W3 one-loop self-energy contributions). Incidentally, our 
parameter choice when illustrating the unconstrained MSSM case below tends to maximize 
the p contributions (due to the simplifying assumption ttiq l ~ ^t R i corresponding to 
m t A t ~ rriQ L ), but this is not a generic situation in the unconstrained MSSM, and can 
therefore be considered as conservative concerning the limitation of present p-parameter 
constraints on such large A t scenarii. Indeed, it is interesting to note that in the minimal 
SUGRA case, additional constraints on the soft-SUSY breaking scalar masses and trilinear 
couplings quite generically imply small Ap(t, b) values. For the typical set of mSUGRA 
parameters chosen in the next sections, the corresponding contributions to the p parameter 



S 



are well below the present exclusion limit, often by an order of magnitude. This can be 
traced back to the renormalization group (RG) evolution [20 1 starting with mSUGRA 
universality conditions at the GUT scale, which indicates that in most cases 3> m t A t 
at the electroweak scale, even when A t is relatively large [] so that we are in the decoupling 
regime for Ap(t,b). 



3. Associated production at proton colliders 



In this section, we examine the associated production of Higgs bosons with a squark pair 



at the LHC. This is an extension of the analysis presented in Ref. |L4| for the associated 
h boson production with stop pairs, with a more detailed discussion of the calculation 
and physics results. In particular we shall illustrate both the case of the unconstrained 
MSSM as well as what is expected in the more constrained minimal SUGRA case. 



3.1 Analytical results 

At lowest order, i.e. at 0{Gpo^), the final state qiq^ where $ is a CP-even Higgs boson 
in pp (or pp) collisions is initiated by the Feynman diagrams shown in Fig. 1. 
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+ 



9 



z _ 



+ 



9 



+cros. 




1 This is actually correlated to the fact that RG evolution indicates that toq l > m tR (or even mg L 3> 
m tR ) in most mSUGRA cases, and also that A t itself cannot be very large, as we shall discuss in more 
details in section 3.2.2 below. 
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Figure 1: Feynman diagrams for the production of a CP-even Higgs boson $ in association 
with a pair of squarks via gg fusion and qq annihilation. 

There are 10 diagrams for the gg fusion mechanism [including those with the quartic 
gluon-squark interaction and the three-gluon vertex] once the various possibilities for 
emitting the Higgs boson from the squark lines and the crossing of the two gluons are 
added, and 2 diagrams for the qq annihilation process^. The 0(Gp) contribution from 7 
and Z-boson exchange diagrams, as well as the flavor changing contribution with a gluino 
exchange diagram are negligible since they are suppressed, respectively, by the additional 
weak coupling factor and the very small mixing between light quarks and top/bottom 
squarks through the gluino interaction, the mixing being due to weak interactions again. 

Due to CP-invariance which forbids the q~iq~iA couplings, the pseudoscalar A boson 
can only be produced in association with q\q\ pairs; this process [which therefore should 
be suppressed by phase space, relative to the production in association with q\q\ pairs] 
has been discussed recently in Ref. |I6[ . 



Due to the larger gluon luminosity at high energies, resulting from the small values of 
the parton momentum fraction x that can be reached and the large probability density 
of the gluon at small x, the contribution of the gg-fusion diagrams is much larger than 
the contribution of the qq annihilation diagrams at LHC energies: at least for most of 
the squark and Higgs boson masses considered in our analysis, it represents roughly 95% 
of the total contribution. In contrast, it is just the reverse at Tevatron energies, where 
the qq initiated process dominates and gives about 80 % of the total contribution. Apart 
form the reduced gluon luminosity at lower energies, this is of course also due to the fact 
that Tevatron being a pp collider, the qq hard process with valence quarks (antiquark) is 
enhanced. 

In terms of the momenta (71,(72 of the initial gluons or quarks, and the momenta 
Pi,P2,P3 of the final two squarks and Higgs boson respectively, the amplitude squared 



2 Actually, in a standard R x type gauge, it is necessary to also add two graphs, corresponding to the 
diagrams of Fig. 1 which involve the triple gluon vertex, but with the initial gluon lines replaced by 
ghosts. Alternatively, one can perform the calculation without ghosts, by choosing e.g. the axial gauge; 
as a cross-check of our calculation we have done it in both gauges and found complete agreement. 
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I-A^lg of the gluon initiated subprocesses gg — > qiq^ [averaged over color and spin], 



calculated with the help of Mathematica , is given by: 

1 



\M\ 



2 
99 



16 



x 



x (Ay/2 G F Mi oL s gt) x — 

64 x 4 i '* i ?i<?i g 

Ci + C 2 Ax + C 3 A 2 + 04^x^2 + C 5 A 2 + C 6 A 2 



(19) 



where g s and Gp are the strong and Fermi coupling constants, and the coefficients Cj are 
given by: 



1 



Cl 

c 2 



(pi.p 2 s - 2x) 



- (pi.p 2 s - 1xfy\V2 + 4(m|s - 4yi)(m|.s - 4y 2 )(:r 2 - 2j/i2/ 2 ) 
2m|s(8pi.g_(p 2 .gi) 2 - 8pi.g!((p 2 .g + ) 2 - 2y 2 ) + p 1 .q + y 2 ) 



+y 2 (-46pi.g_(p 2 .gi) 2 + 7sx + 46pi.gi((p 2 .g+) 2 - 2y 2 ) + 10pi.g+j/ 2 ) 



c 3 



(Pi-P2« ~ 2a;) 
4S 2 y 2 y 2 



2m?s(8p 2 .g_(pi.gi) 2 - 8p 2 .gi((pi.g+) 2 - 2j/i) + p 2 .q+y 1 ) 



+y 1 (4:6p 2 .q^(p 1 .q 1 ) 2 + 7sa; - 46p 2 .qi((pi.q + ) 2 - 2y x ) - 10p 2 .q + y 1 ) 



{pi.p 2 s -2x) 2 7 



4 m 



7 



s 2 ym 



+ t;Pi-1+(-- + — 1 ) + ^h— (9pi.q-P2-q- ~ ^P2-q+s){m q s - Ay 2 
ys z 2 s y\ s z y 2 2 

1 

{{y - 8x)(pi.p 2 s - 2x) 2 + 9m|pi.g_p 2 .<?_% 2 ) + 7 



c 5 



2S 2 y 2 



2m 2 ~i 2 (4(p 2 .g + ) 2 - 9y 2 ) + (92(p 2 .g + ) 2 + 28p 2 .q + s + 7s 2 - UAy 2 )y 2 



+mq.sy 2 (-46(p 2 .q+) - 7p 2 .q + s + 72y 2 ) 



c 6 



2s 2 y\ 



2mls 2 (A( Pl .q + y - 9 Vl ) + (92( Pl .q + ) 2 - 28 Pl .q + s + 7s 2 - UAy^yf 



+m 2 =sy 1 {-46(pi.q + ) 2 + 7p x .q + s + 72y x 



(20) 



In the previous equations, with the help of q± = q± ± q 2 , we have used the following 
abbreviations to simplify the result: 



At 



1 



[s + 2p 2 .q + ) 



A, 



1 



(s - 2p x .q^ 



and 



x = pi.q 2 p 2 .qi + pi.qip 2 .q 2 , y { = Pi.q1Pi.q2 , y = P1.q2P2.q2 



(21) 
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In the case of the amplitude squared |.M|L of the quark initiated subprocesses qq — > qi.qi<& 
[averaged over color and spin] , the analytical expression is much simpler then the previous 
expression for |.M| 2 and is given by: 

M« = 9^4 ^ G FM*glMd) * 2 >< J2^ + A l d * + ^A 1 A 2 d a ) (22) 

where A x and A 2 are given previously, and in terms of Ri = 2p 2 + q + and R 2 = 2p 1 — q + , 
the coefficients dj are given by 

d x = 2(R 1 ■ q^ ■ q 2 ) - S -R\ , d 2 = 2(R 2 ■ qi )(R 2 ■ q 2 ) - S -R\ , 

ds = (Ri ■ qi){Rz ■ ft) + (^2 ■ qx){Ri ■ ft) - • R 2 ) (23) 

The hard scattering process is convoluted with the initial parton densities, and inte- 
grated over the final state phase space. The resulting q~iq~i§ cross-section reads 

I dT dx r 
o{pP -> <Mi$) = — / — — dPS(q 1 +q 2 ;p 1 ,p 2 ,p 3 ) x 

S^2[q(x,Q f )q(r/x,Q f ) + q^q\\M\ 2 q , + p(x, Q/)<?(t/z, Q/)|.M|1,} (24) 

where q(x,Qf), q(x,Qf) [running over the five flavors w, d, c, s, 6] and g(x,Qf) designate 
respectively the quark, antiquark and gluon density functions with momentum fraction x 
and factorization scale Qf. In eq. (|24]), s is the hadron cm. energy squared and as usual 
s = x\x 2 s = ts is the parton cm. energy, with r min = (2m^ + M$) 2 /s. In eq. (|24D (iPS 
(<Zi + Q.2\PiiP2iPz) is an element of the 3-body final state phase space normalized as 



JdPS = J(27i) i 5\q 1 + q 2 -p 1 -p 2 -p 3 ) 



d 3 pi d 3 p 2 d 3 p 3 



(2tv) 3 2E 1 (2tt) 3 2£ 2 (2tt) 3 2E 3 
d(cos6) d<f> dE x \ dE l2 , (25) 



(27I") 4 8 J-l JO J_Ei,min jE 12 ,mm 

where E\, E 2 are the energies of the two produced squarks with E\ 2 = (E\ — E 2 )/2, 
and £"3 the energy of the Higgs particle $. The remaining integral in eq. (p5j) is over 
appropriately defined angles 9 and 0, describing the motion with respect to the beam axis 
of the 3-momenta pi, p 2 of the two produced squarks^. The integration bounds are given 
from kinematics as follows: 

s + m| - (m§ 4 + M$) 2 

-El,min = Tflqi , -El, max = " 7^ (26) 

2y s 



3 Actually four angles would be needed, but one angular integration is eliminated from energy conser- 
vation and another (azimuthal) angle integration simply gives a factor of 2ir included in eq. (|25|). 
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whereas the expressions for i?i2,mm and -Ei2,max ar e more involved: 

-b =F \/b 2 - Aac 



- C/ 12,min/max 

with a, b, and c given by: 

a = 2EiVI — s — nig. , I 



2a 



(27) 



£1) (M, 



2 2 
Hi 



V§ - Ei) -4ml 



El + Ml- 2ml 



(28) 



Finally, the resulting six-dimensional integral eq. (|24|) over the remaining phase space and 
over the parton luminosities is performed numerically with the standard Vegas Monte- 
Carlo integration routine [28|. 



3.2 Numerical illustrations 

We now turn to the numerical results, that we will use to illustrate the production of 
the lightest h boson in the decoupling limit [with Ma — 1 TeV] in association with the 
lightest top squarks. Throughout the analysis we will use the most recent CTEQ4 pa- 
rameterizations of the parton density functions [29|| , at leading order, with a factorization 
scale Qf = a/I. Note that we also use a (one-loop) running a s , defined at the scale s. 
The top quark mass is fixed to m t = 175 GeV. 



3.2.1 Unconstrained MSSM 

To simplify the phenomenological analysis in the case of the unconstrained MSSM, we will 



m, 



assume the left- and right-handed stop mass parameters to be equal, m^ L = m^ R 
For illustration, we have chosen two somewhat extreme values for tan/5: tan/5 = 3 and 
tan j3 = 30. Large values of tan (5 maximize the h boson mass and therefore circumvent 



the present experimental lower bound on in the decoupling limit, Mh > 90 GeV [IT], 
almost independently of the values of the parameter A t chosen in our illustrations. How- 
ever, large tan (3 values tend to worsen the constraints from CCB and the p parameters, 
as has been discussed in previous subsection 2.2. For smaller values of tan/3, the 
coupling is simply obtained by rescaling the parameter A t = A t — \ij tan (3. An advantage 
of low tan (3 values is that one can then reach relatively large values of A t and thus smaller 
m~ tl before encountering problems with the CCB constraints and/or the rho parameter. 
Mh also becomes smaller, allowing for more phase space and thus increasing the cross 
section, but eventually reaching its present lowest limit for a too small value of A t . 

In Fig. 2, the pp — > t]t\h cross section [in pb] is displayed as a function of the lightest 
t mass for the value tan (3 = 30 in the case of no mixing A t ~ [A t = 0, /i = —100 GeV] 
and moderate mixing [A t = 500 GeV and /x = —100 GeV]; and for the value tan/3 = 3 in 
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the large mixing case A t ~ 1.4 TeV [A| = 1.2 TeV and // = —600 GeV]Q. For comparison 
the cross section for the standard-like pp — > tt/i process [that we indeed recalculated 
independently as a cross-check of our numerical procedure] is of the order of 0.5 pb for a 
Higgs boson mass Mh — 100 GeV | 17fl . The cross section behaves as follows: 



a(pp -!• tJih) [pb] 
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Figure 2: The production cross section a{pp — > titih) [in pb] at the LHC as a function of 
the ti mass and three different choices of the other parameters: tan/3 = 30, A t = (0.5) 
TeV; tan/? = 3, A t = 1.2 TeV. 



- In the case where there is no mixing in the stop sector, t\ and t 2 have almost the 
same mass [which, up to the small contribution of the D-terms, is constrained to be larger 
than m^ s ~ 165 GeV for m~ > 0] and approximately the same couplings to the h boson 
since the m^/Mf components in eq. (12) are dominant. The cross section in Fig. 2, which 
should then be multiplied by a factor of two to take into account the production of both 
top squarks, is comparable to the cross section for the SM-like process pp — > tth in the 
low mass ranged mi < 200 GeV. 

- For intermediate values of A t the two components of the htit\ coupling interfere 

4 Note that for given A t and m; values, the cross sections in Fig. 2 are slightly smaller than the 
corresponding ones given in Ref. |14| , where Mh was a little lower than the present value due to a 
different choice of the parameter tan (3. 

5 In scenarii where the t masses are related to the masses of the light quark partners rriq, the mass 
range mi < 200 GeV in the no mixing case is, however, ruled out by present experimental constraints 
on rriq [which approximately corresponds to the common mass of the squarks of the first two generations 
and the sbottom] at the Tevatron [O . 
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destructively [with our conventions, sin2# t < in the relevant parameter space, see eq. 
(f|)] partly canceling each other and resulting in a small cross section, unless < 100 
GeV. For some value of At } 9htih ~ ^ and the cross section vanishes. 

- In the large mixing case, A t ~ 1.2 TeV, a(pp — > tit\h) can be quite large without 
conflicting with the present bounds on M^, the rho parameter, or the CCB constraints 0. It 
is above the rate for the standard process pp — > ith for values of mt smaller than 210-220 
GeV, approximately. If t\ is lighter than the top quark, the t\t\h cross section significantly 
exceeds the one for ith final states. For m f - = 150 GeV for instance, a(pp — > t\t\h) is 
about a factor five larger than a(pp — > tih). 



a(pp -» Utth) [pb] 




100 300 500 700 900 1100 1300 

At (GeV) 

Figure 3: The cross section a(pp — > t-\t\h) [in pb] at the LHC as a function of /l t for fixed 
m ?l =165 GeV, tan/? = 3; and 250 GeV, tan/3 = 30. 

In Fig. 3, we fix the lightest top squark mass to mj = 165 and 250 GeV, and display 
the pp — > gg + qq — ► titi/i cross section as a function of A t . The value of tan/3 and 
are fixed to tan/3 = 3, /i = —600 GeV and tan/3 = 30, /i = —100 GeV respectively, in 
order to illustrate that the A t range allowed from the different constraints also depends 
somewhat on the values of tan/3, as already explained above. For = 165 GeV, the 
production cross section is rather small for the no mixing case and even smaller for the 
intermediate mixing [becoming negligible for A t values between 50 and 150 GeV], and 
then becomes very large, exceeding the reference cross section [a(pp — > ith) ~ 0.5 pb for 
M/i ~ 100 GeV] for values of A t above ~ 0.9 TeV. For fixed value of A t , the cross section 

6 The range of i% masses below < 130 GeV where the cross section is the largest, gives too large 
contributions to the p parameter for this specific large value of At, and is also excluded already by the h 
boson mass lower bound. 
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decreases with increasing t\ mass since the phase space becomes smaller. It is however 
still sizeable for m k = 250 GeV and A t > 1.4 TeV. 

Note that for fixed t mass and coupling, the cross section becomes smaller for larger 
values of A t , if A t < VGrriq, because Mh increases [|J and the process is less favored by 
phase-space; in the reverse situation, A t > VGrriq, the h boson mass will start decreasing 
with increasing A t [reaching values below Mh < 90 GeV when e.g. A t ~ 1.3 TeV for 
mi = 165 GeV and tan/3 = 3, where we thus stopped the corresponding plot] and the 
phase-space is more favorable to the reaction. 



3.2.2 The mSUGRA case 



Let us now turn to the special case of the minimal SUGRA model [Tj|, in which universal- 
ity at the GUT scale implies that the only free parameters [l[] are the values of mo, 771-1/2, 
Aq, respectively the common scalar mass, gaugino mass and trilinear scalar coupling, plus 
the sign of the /i parameter and tan (3. All the physical parameters, and in particular the 
squark masses, are obtained from a renormalization group (RG) running from the GUT 



scale to the weak scale [20], where the electroweak symmetry breaking constraints fix the 



\n\ and B parameters of the Higgs potential at the electroweak scaleQ. 

For our numerical illustrations, Fig. 4, we choose specific values of the above param- 
eters in such a way that one of the top squarks can be sufficiently light in order to have 
sizeable cross sections. In general, the cross section in the mSUGRA scenario may be as 
substantial as in the unconstrained MSSM cases illustrated above, but this occurs in a 
more restricted region of the parameter space. This is essentially due to the fact that it 
is generically very difficult to have almost degenerate ti and tn in mSUGRA as already 
discussed in section 2 and as is clear from eqs. (^|), so that the stop mixing angle which is 
controlled by the ratio A t /{rhj — m~ ) can become large only for very large A t . Moreover 
considering the RG equation driving the A t parameter [|20J : 



A A 1 3 

= /3(A) = g^(E C ^M, + QY t 2 A t + Y b 2 A b ) (29) 

where t = InQ, g iy Yi are the gauge and Yukawa couplings and M, the gaugino masses 
[with the constants q given by: c\ = 26/15, c 2 = 6 and c 3 = 32/3], one can see that for 
most cases | A t | tends to decrease when the energy scale Q is decreasing from GUT to low- 
energy. Indeed, the q's in eq. (p9|) are all positive and in mSUGRA Mj(GUT) = my 2 > 0, 
so that 13 {A t ) > for A t > 0. This accordingly makes a large A t value at low energy less 
likely, since Aq = A t {GUT) would have to be even larger, which may conflict with the 
CCB constraints among other things, as it was discussed in the previous section. The 

7 The RG evolution of the relevant parameters and calculation of the physical masses and couplings is 
done with a numerical code |3C}| including in particular supersymmetric threshold effects at the one-loop 
level, and a consistent treatment of the electroweak symmetry breaking. 
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only way to have an increasing | A t | when running down to low energy is if A < with A Q 
small enough that (3{A t ) in eq. ( ^9|) remains positive, but the latter positiveness requires 
a relatively large value, which in turn implies from eq. (|6]) that m t - cannot be small. 
[Note that, due to the universality of the trilinear scalar coupling, the above features are 
qualitatively unchanged even when the Yukawa bottom contribution in eq. ( p9|) is not 
neglected.] 




Figure 4: The cross section a(pp — > titih) [in pb] at the LHC as a function of the t\ mass 
in the mSUGRA case: tan/5 = 2.5(30), m 1/2 = 300 GeV, A = -2 TeV, sign(//) = +. 

The mSUGRA cross section as illustrated in Fig. 4 actually corresponds to a situation 
where \A t \ at the relevant electroweak scale is not very large, \A t \ < 1 TeV, while the 
effective coupling \A t \ becomes larger, at least for the small tan/? =2.5 case, because /i 
[which we determine consistently from radiative electroweak symmetry breaking] turns 
out to be rather large, fi ~ 1 TeV. In this case, for a sufficiently low m~ tl value the 
cross section is potentially as large as in the case of the unconstrained MSSM, but this 
situation clearly occurs in a restricted region of the mSUGRA parameter space. On the 
other hand, as discussed in subsection 2.2, the mSUGRA case is accordingly less subject 
to the CCB constraint, since A t is not too large, and also no restrictions are imposed by 
the p parameter until the t\ mass becomes really small. However, one should take into 
account another constraint, namely that of requiring the lightest neutralino [rather than 
the lightest top squark] to be the lightest SUSY particle (LSP). This is the reason why 
in Fig. 4 the plots start at m~ tl values larger than m x o ~ 120 GeV, which corresponds to 
our choice of the basic mSUGRA parameters. 
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3.2.3 Associated production at the upgraded Tevatron 

For completeness, we illustrate some values of the t\t\h production cross section at the 
upgraded Tevatron with y/s = 2 TeV. Naively, one may expect that such a three particle 
final state process will be very disfavored by phase space in the parton frame, except for 
a very light t± and very light Higgs bosons which are already excluded by present experi- 
mental constraints. However, this clear phase space suppression is partially compensated 
by the fact that the qq hard scattering process [which scales as the inverse of the parton 
cm. energy s] dominates the cross section at the Tevatron. More precisely, due to the fact 
that Tevatron is a pp collider, the cross section benefits form the larger quark/antiquark 
densities relative to the LHC case. 




Figure 5: The production cross section a(pp — > tit\h) at the Tevatron [in pb] as a function 
of the t\ mass for A t = 1.2 TeV, /i = —600 GeV and tan/3 = 3. 

The cross section a{pp — > ii^h) is shown in Fig. 5 at a cm. energy */s = 2 TeV 
in the case of large A t values, since otherwise the cross section is clearly too small for 
the process to be relevant. For the planned luminosities of / Cdt ~ 10 fb _1 ||, one can 
see from Fig. 5 that up to a few tens of events might be collected for m~ tl ~ 130 GeV 
[smaller masses for such a large A t value, are already ruled out by the present constraint, 
as discussed above], before efficiency cuts are applied. This gives some hope that this 
process might be useful at the Tevatron too. We shall discuss in the next subsection the 
possible signal for this process, which is similar for the LHC and Tevatron. 
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3.2.4 Signal 

Finally, let us discuss the signal for the pp — ► t\t\h process. If the mass of the lightest 
top squark is only slightly larger than that of the lightest neutralino [which is expected 



to be the LSP], the dominant t\ decay channel [31| will be the loop induced decay into a 
charm quark plus a neutralino, t\ — > c%°- F° r larger mi x values, the lightest top squark 
can decay into a b quark and a chargino, t\ — > b\ + and possibly into a t quark and the 
LSP, t — > txi- [We assume that the strong decay into gluinos does not occur]. In the 
interesting region where the cross section a(pp — > t±tih) is large, i.e. for relatively light t±, 
the decay mode t% — > b\ + can be dominant, unless the mass difference — m x + is very 
small. Assuming that the partners of the leptons are heavier than the lightest chargino, 
Xi w iH mainly decay into the LSP and a real or virtual W boson, leading to the final 
state 

t\ — ► bW + + missing energy (30) 

This is the same topology as in the case of the top quark decay, t — > bW + , except that in 
the case of the top squark there is a large amount of missing energy due to the undetected 
LSP. If sleptons are also relatively light, charginos decay will also lead to lv\\ final states. 
The only difference between the final states generated by the tth and tih processes, will 
be due to the softer energy spectrum of the charged leptons coming from the chargino 
decay in the former case, because of the energy carried by the invisible neutralinos. 

The Higgs boson can be tagged^ through its h — > 77 decay mode . In the decoupling 
limit, and for light top squarks and large A t values, the branching ratio for this mode can 
be enhanced [up to ~ 20%] compared to the 77 branching ratio of the SM Higgs boson 



13], because of the additional contributions of the t-loops which interfere constructively 



with the dominant H^-loop contribution. 

Although a detailed analysis taking into account backgrounds and detection efficien- 
cies, which is beyond the scope of this paper, will be required to assess the importance 
of this signal, it is clear that in some areas of the MSSM parameter space, 77+ charged 
lepton events can be much more copious than in the SM, and the contributions of the 
pp — > tth process to these events can render the detection of the h boson much easier than 
with the process pp — > tih alone. This excess of events could then allow to measure the 
htiti coupling, and open a window to probe the soft-SUSY breaking parameters of the 
stop sector. In particular one can obtain some information on the trilinear coupling A t , 
which is rather difficult to extract from the study of the production of SUSY particles at 
the LHC as shown in Ref. |33| . 



^Decays of the h boson, produced in association with tt pairs, into the main decay channel h — > bb 



have also been discussed in the literature 1 32 
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4. Associated production at e + e colliders 



At future linear e + e~ colliders, the final state qiq±h may be generated in three ways: (i) 
two-body production of a mixed pair of squarks and the decay of the heaviest squark 
to the lightest one and a Higgs boson [Fig. 6a]; (ii) the continuum production in e + e~ 
annihilation^ [Fig. 6b]; and (Hi) the continuum production in 77 collisions [Fig. 6c]. The 
same picture holds for the associated production of the heavier CP-even Higgs boson H 
with q~iq~i pairs. In the case of the pseudoscalar boson A only the processes (i) and (ii) are 
possible: because of CP-invariance, the coupling Aq~iq~i is absent and the final state q\q\A 
cannot be generated in 77 collisions at tree-level; furthermore, in e + e~ annihilation the 
graph where the A boson is emitted from the Z line is absent, but there is another graph 
involving the ZAh(H) coupling [see Fig. 6b]. For the associated production of CP-even 
Higgs bosons with g 2 g 2 states, only the processes (ii) and (Hi) are relevant. 



9 Recently, the continuum production in e + e annihilation has also been evaluated in Ref. |l5| using 
the numerical graph calculation program GRACE-SUSY ]34|j , 
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Figure 6: Feynman diagrams for the production of a Higgs boson $ in association with a 
pair of squarks in e + e~ and 77 collisions. 
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4.1 Two— body production and decay 

If the center of mass energy of the e + e~ collider is high enough, one can generate a qiqih 
final state by first producing a mixed pair of squarks, e + e~ — > q~iq~2, through the exchange 
of a virtual Z-boson, and then let the heaviest squark decay into the lightest one and 
the Higgs boson, q 2 —>■ q±h, if the splitting between the two squarks is larger than M^; 
Fig. 6a. For H and A bosons final states, a larger splitting in the stop sector is required. 

The total cross section for the production of two squarks with different masses in e + e~ 
annihilation, e + e~ — > q~iq 2 , including the factor of two to take into account the charge 
conjugated states is given by fl35 



a(e + e -> gig 2 ) = - (a e + f e ) o 12 _ ^ / (31) 

with s the center of mass energy of the e + e~ collider; o"o = 47ra: 2 /3s the QED point cross 
section and A the usual two-body phase-space function, X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + 
xz + yz). The Zq\q\ coupling a 12 is proportional to sin 26^, eq. (15), and in the case where 
the mixing between squarks is absent, this coupling is zero and the cross section vanishes. 
[The QCD corrections to this process, including the SUSY-QCD contributions with the 



exchange of gluinos, are available [35] and can be included.] 

The partial decay width of the heaviest squark into the lightest squark and a neutral 
Higgs boson $ = h, H, A is given by |35 



Tfe - m = X 1/2 (ml, Ml ml) g%- qih (32) 



<I2 



where the two-body phase space function A has been given previously. [Here also the 
QCD corrections are available |J7| and can be included]. As can be seen from inspection 
of eq. fll2D , in contrast with the g^q^ coupling, Ig*^! cannot be too large in most of 
the parameter space, and in fact this coupling is well dominated by its second term, 
proportional to cos2#g, except for very large A t and/or when the mixing between squarks 
is maximal, \6 q \ = 7r/4, in which case it is anyway practically vanishing. 

The branching ratio for the q 2 — > qi<& decay mode is obtained by dividing the partial 
width by the total q 2 decay width. The latter is obtained by summing the widths of all 
possible decay channels. For instance, in the case of the more phase-space favored decay 
channel t 2 — ► tih, one has to include: decays of the t 2 into a bottom squark and a charged 
Higgs or W boson, decays into t% and a neutral Higgs or Z boson, decays into a b quark 
and charginos, and decays into a top quark and neutralinos [or gluinos]. 

In principle, if phase-space allowed, the cross section for the two-body production 
process times the branching ratio for the two-body decay, should be large enough for 
the final state to be copiously produced. However, as discussed above, in a large part of 
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the parameter space the cross section times branching ratio will be roughly proportional 
to sin2#g x cos2^ g , thus giving small production rates in the no mixing and maximal 
mixing scenarii. In addition, the decay width q 2 — > hq\ is in general much smaller than 
the q 2 decay widths into chargino and neutralinos, leading to a small branching ratio. 
Nevertheless, there are regions of the MSSM parameter space where the combination 
sin2# g x cos2#g can be maximal, which occurs typically for a not too small m^-m^ 
splitting and a moderate A t , as can be seen from eq. @. In this case [which is often 
realized in particular in the mSUGRA case, as we shall discuss next] the resonant t 2 
process may dominate over the non-resonant qqh production and the corresponding rate 
is visible for the high luminosities / £dt ~ 500 fb -1 expected at linear colliders [BH]. 

Such a situation is illustrated in the case of the t±tih final state in Fig. 7, where the 
cross section e + e~ — > tit 2 times the branching ratio BR(i 2 — ► t\h) is shown as a function 
of the t\ mass at a cm. energy of \fs = 800 GeV [full lines]. We have chosen a mSUGRA 
scenario with tan/3 = 30, m-i/ 2 = 100 GeV, A = —600 GeV and sign(/x) = +. [The 
broken lines show the contribution of the non-resonant contributions in this case, which 
will be discussed in more detail later]. As can be seen, the cross section can reach the 
level of 1 fb for relatively small m~ tl values, leading to more than one thousand events in 
the course of a few years, with the expected integrated luminosity of / Cdt ~ 500 fb -1 
Moreover, for the same reasons that have been discussed in previous pp section for 



the mSUGRA case, there are no other constraints e.g. from CCB or the rho parameter, 
for the values chosen in our illustration. 
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Figure 7: The production cross section cx(e + e — > t^ih) [in fb] as a function of in the 
mSUGRA case; tan/5 = 30, m 1/2 = 100 GeV and A = -600 GeV. 
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As discussed in the previous section the top squarks in this mass range will mainly 
decay into a charm+neutralino, t\ — > cx±, or a b quark and a chargino, ti — > b\ + ■ In 
this latter case the lightest chargino, xt will decay into the LSP and a real or virtual W 
boson, leading to the same topology as in the case of the top quark decay, but with a 
large amount of missing energy due to the undetected LSP. However, in e + e~ collisions 
one can use the dominant decay mode of the lightest Higgs boson, h — > bb. The final state 
topology will consist then of 46 quarks [two of them peaking at an invariant mass M^, 
which would be already measured by means of another process], two [real or virtual] W's 
and missing energy. With the help of efficient micro-vertex detectors, this spectacular 
final state should be rather easy to detect in the clean environment of e + e~ colliders. 



4.2 Production in the continuum in e + e collisions 

There are three types of Feynman diagrams leading to final states with qiqi pairs and a 
CP-even Higgs boson h or H, in e + e~ annihilation: Higgs boson emission from the q± 
states [which are produced through s-channel photon and Z-boson exchange] , Higgs boson 
emission from the q 2 state [which is produced with q~\ through Z-boson exchange], and 
Higgs boson emission from a virtual Z boson which then splits into q\q\ pairs. When the 
Higgs-gigi couplings are large, the two last types of Feynman diagrams give negligible 
contributions since the virtuality of g 2 is large and the Zq\q\ coupling not enhanced. 
[In the opposite situation where the squark mixing angle 6 q is non-zero but moderate, 
when the (resonant) q 2 contribution largely dominates the cross section, it is then very 
well approximated by the on-shell q\q 2 production with subsequent decay q\ — * Qih, as 
illustrated in the previous subsection 4.1 in the mSUGRA case]. For the CP-odd Higgs 
boson the first and last possibilities are absent because of CP-invariance which forbids 
the couplings of A to gauge bosons and qiq~i pairs. 

Let us first discuss the process e + e~ — > qiq^ with <3> the CP-even Higgs boson h or 
H, and q\ any of the two squarks. In terms of the scaled energies of the two final state 
squarks x lj2 = 2E 1:2 /y/s, the Dalitz plot density for the reaction reads 



da 



e e 



qiqi®) = (7 



.2 
W 



(33) 



fc=l,6 



Using the reduced masses //j = m~./s,f/,j = m~ /s [qj is the virtually exchanged partner 
of the produced squark qi\, //$ = M|/s and z = M^/s, and the couplings given in section 
2, the various amplitudes squared are as follows: 
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In these equations, 2/1,2 = 1 —2:1,2, = 1 — 2/1 — 2/2 + Ai* — 2/!« and the (reduced) couplings 
g<s> of the $ boson to two Z bosons [in units of eMzf(swCw)] are given by: = sin(/3 — a) 
and = cos(/3 — a). 

For the final state $qiqi in the case where the g^q^ coupling is large [leading to 
a large splitting between the q~\ and q\ states], the virtuality of the squark q\ is large 
and the Zq\q\ coupling not enhanced; the Dalitz plot density of the process can be then 
approximated by the very simple form 
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For the associated production of the pseudoscalar Higgs boson with a pair of squarks, 
e + e~ — > q~iq~iA, the Dalitz plot density has a much simpler form since only the diagrams 
with the virtual q~j exchange and the one involving the ZAh(H) vertices, are present. It 
reads: 



da 

dxix 2 



[e + e -> qiqiA) = a 



(a 2 e+V 2 e) 



(36) 



fc=l,3 



where the amplitudes squared, using the same notation as previously, are given by 
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with the Z&A couplings gzAh = cos(/3 — a)/ (2swCw) and Qzah = — sin(/3 — a)/ (2swCw) 
(in units of e); here again qj is the virtually exchanged partner of the produced squark & 
with qj j£ q~i. 

To obtain the total production cross section a(e + e~ — ► q~iq~i§), one has to integrate 
over the scaled variables = 2E\^j \fs in a way totally similar to the one described in 
previous section 3, see eqs. (p5|-[28|). 



a(e + e- -> titjh) [fb] 




0.01 

280 300 
mi, (GeV) 

Figure 8: The cross section <r(e + e~ — > t\t-Ji) [in fb] as a function of the ti mass and three 
different choices of the other parameters: tan/3 = 30, A 4 = (0.5) TeV; tan/3 = 3, A t = 
1.2 TeV. 



To illustrate the magnitude of the cross sections, we show in Figs. 8-9 the rates for the 
tit\h final state in the decoupling limit at y/s = 800 GeV and tan/3 = 30(3) as a function 
of the ti mass for A t ~ and A t = 0.5 TeV (A t = 1.2 TeV) respectively [Fig. 8], and as 
a function of A t for m h = 165 GeV (tan/3 = 3), and 250 GeV (tan/3 = 30) [Fig. 9], i.e. 
for the same scenarii as in Figs. 2 and 3 of section 3 [the discussions on the constraints 
on the parameters A t and m~ tl given there apply also in this case]. 

As can be seen the trend is similar to the one discussed in the previous section for the 
associated production at proton colliders. For not too large t\ masses and large values of 
the parameter A t , the production cross sections can exceed the value <r(e + e~ — ► t\tih) ~ 1 
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fb. [Note that the cross section for the SM-like process e + e — > tth |18j is of the order 



of 2 fb for Mh ~ 130 GeV at a cm. energy ^fs = 800 GeV.] This provides more than 
one thousand events in a few years, with a luminosity / Cdt ~ 500 fb -1 : a sample which 
should be sufficient to isolate the final state [with the topologies discussed in section 4.1] 
and measure the gi x i xh coupling with some accuracy. 

At lower energies, the production cross section will be as large as in the previous case 
for small t\ masses, but will be limited by phase space for larger values. This is illustrated 
in Fig. 10, with the same choice of parameters as in Fig. 8, but with a cm. energy of 
sfs = 500 GeV. Thus for very large values of A t , the cross section cr(e + e~ — > tii\h) can 
exceed the fb level for stop masses below < 150 GeV even at a 500 GeV collider. With 
~ 500 fb -1 luminosity, this would hopefully allow to detect the final state and measure 
the coupling g^f j with some accuracy. In the mSUGRA case the cross section generally 
follows also the same lines as in section 3.2.2 and Fig. 4, i.e. it can be as large as in the 
case of the unconstrained MSSM, but in relatively smaller area of the SUSY parameter 
space, for reasons that have been given earlier. 

Note finally that in the unconstrained MSSM case, the continuum production cross 
section in e + e~ annihilation cr(e + e~ — > t\tih) is often larger than the resonant cross 
section for the production of tit 2 and the subsequent 2-body decay t 2 — > t\h, but this is 
not generic. Indeed, in a situation where both a non-negligible m^-m^ splitting and a 
moderate A t occurs, provided there is sufficient phase space allowed, the production via a 
resonant t 2 becomes competitive and even dominant, as illustrated in a typical mSUGRA 
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o(e+e- -» hhh) [fb] 




m fl (GeV) 

Figure 10: Same as Fig. 8 but with y/s = 500 GeV. 
case in the previous subsection 4.1. 

4.3. Production in the continuum in 77 collisions 

Future high-energy e + e" linear colliders can be turned into high-energy 77 colliders, with 
the high energy photons coming from Compton back-scattering of laser beams. The cm. 
energy of the 77 collider is expected to be as much as ~ 80% of the one of the original 
e + e~ machine. However, the total luminosity is expected to be somewhat smaller than the 
one of the e + e~ mode, leading to a smaller number of events for the same cross section. 

The final state q^ih or H can be generated by emitting the h, H bosons from the 
squark lines in the process 77 — > 5^; Fig. 6c. The pseudoscalar Higgs boson A cannot 
be produced in two-photon collisions in association with a pair of squarks because of 
CP-invariance; the production of the A boson is only possible in association with qiq 2 
states which is less favored by phase space. 

The differential cross section of the subprocess 77 — > with $ a CP-even Higgs 

boson is given by 

Ar( 77 - M&) = Ys XdFSx |M| ^ (38) 

where dPS is the same element of the 3-body phase space discussed in section 3 for 
hadron colliders, and |M|L, is the amplitude squared of the subprocess; taking the same 
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convention as in section 3 for the momenta of the initial and final particles, and averaging 
over spin, the latter is given by 



|M| 
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77 



- e\ (AV2 G f M| g% Mi e 4 ) \ 9l + g 2 A x + g 3 A 2 + g 4 A 1 A 2 + g h A\ + g 6 A*] (39) 



where A lj2 have been defined in eq. (19); in terms of the variables yi t2 , x i V defined also in 
eq. (19), and s, the cm. energy of the 77 subprocess and the squark mass respectively, 
the coefficients ^ read: 
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(40) 



Integrating over the phase-space in a similar way as the one discussed in section 3, one 
obtains the total production cross section of this subprocess. One can then convolute with 
the photon spectra [some examples of spectra can be found in Ref. ]3£| for instance] to 
obtain the final cross section 17(77 — > QiQi$)- 111 the following, we will not use any photon 
spectrum for simplicity, we will just exhibit and discuss the production cross section for 
the subprocess. 

The total cross section for the process 77 — > txt\h is shown in Fig. 11 at a two- 
photon cm. energy y^ 77 ~ O-^V^ee = 600 GeV and as a function of the t\ mass, without 
convolution with the photon spectrum and with the same inputs and assumptions as in 
Fig. 8 to compare with the e + e~ mode. Because the cm. energy of the 77 collider is only 
~ 80% of the one of the original e + e~ machine, the process is of course less phase-space 
favored than in the e + e~ mode. Nevertheless, the cross section for the t-J;\h final state is 
of the same order as in the e + e~ mode for cm. energies not too close to the kinematical 
threshold, and the process might be useful to obtain complementary information since it 
does not involve the Z-boson and t 2 exchanges. If the luminosities of the 77 and e + e~ 
colliders are comparable, a large number of events might be collected for small stop masses 
and large A t values. 
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Figure 11: The cross section (7(77 — > titih) in 77 collisions [in fb] at a center of mass 
energy ^/s 77 = 600 GeV as a function of . The other parameters have the same values 
as in Fig. 10. 

5. Conclusions 

We have calculated the cross sections for the production of neutral Higgs particles in as- 
sociation with the supersymmetric scalar partners of the third generation quarks, e + e _ — > 
q~iq~i&, at future high-energy hadron colliders and at e + e~ linear machines, in both the 
e + e~ annihilation and the 77 fusion mode. Complete but rather simple analytical for- 
mulae for the squared amplitudes of these three-body processes are given in the case 
where the two final state squarks have equal mass. In the case of e + e~ collisions, the 
production of both the CP-even and CP-odd Higgs bosons were calculated, while in pp 
and 77 collisions only the production of the CP-even Higgs particles has been discussed, 
since pseudoscalar Higgs bosons cannot be produced at tree-level with a pair of equal 
mass squarks at these colliders. 

In the framework of the Minimal Supersymmetric extension of the Standard Model, 
we have then investigated the magnitude of the cross sections at the upgraded Tevatron 
and the LHC for hadron machines, and at a future e + e~ linear machine with cm. energies 
in the range =500-800 GeV; in the latter case, the 77 option of the collider has also 
been considered. In this numerical analysis, we focussed on the production of the lightest 
CP-even Higgs boson in the decoupling regime, in association with a pair of the lightest 
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top squarks. Final states with the other Higgs bosons or with qiq 2 pairs should be phase- 
space suppressed. 

At the LHC, the production cross sections can be rather substantial, especially in the 
case of rather light top squarks, < 250 GeV and large trilinear coupling A t > 1 TeV. 
In this case the production rates can exceed the one for the standard-like production of 
the h boson in association with top quarks, pp — > tih, which is expected to provide a signal 
of the h boson in the 77/ ± channel. Although a detailed Monte-Carlo analysis, which is 
beyond the scope of this paper, will be required to assess the importance of this signal 
and to optimize the cuts needed not to dilute the contribution of the i\i\h final states, 
it turns out that in a substantial area of the MSSM parameter space, the contribution of 
the top squark to the signal can significantly enhance the potential of the LHC to 
discover the lightest MSSM Higgs boson in this channel. 

At e + e~ colliders with cm. energies around ^/s = 800 GeV and with very high 
luminosities / Cdt ~ 500 fb _1 , the process e + e~ — > tit±h can lead to several hundreds 
of events, since the cross sections [in particular for rather light top squarks, < 250 
GeV and large trilinear coupling, A t > 1 TeV] can exceed the level of a 1 fb, and thus 
the rate for the standard-like process e + e~ — > tih. In the case where the top squark 
decays into a b quark and a real/virtual chargino, the final state topology [with 46 quarks, 
missing energy and additional jets or leptons] will be rather spectacular and should be 
easy to be seen experimentally, thanks to the clean environment of these colliders. In 
the 77 option of the e + e~ collider, the cross sections are similar as previously far from 
the particle thresholds, but are suppressed for larger masses because of the reduced cm. 
energies. For 77 luminosities of the same order as the original e + e~ luminosities, the titih 
final state should also be observable in the two-photon mode, at least in some areas of 
the MSSM parameter space. 

The production cross section of the titih final state is directly proportional to the 
square of the iiiih couplings, the potentially largest electroweak coupling in the MSSM. 
Analyzing these final states at hadron or electron-positron machines will therefore allow 
to measure this coupling, opening thus a window to probe directly some of the soft-SUSY 
breaking scalar potential. 
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